Quantum depinning of the magnetic vortex core in micron-size permalloy disks 
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The vortex state, characterized by an in-plane closed flux domain structure and an out-of-plane 
magnetization at its centre (the vortex core), is one of the magnetic equilibria of thin soft ferromag- 
netic micron-size dots. In the last two decades many groups have been working on the dynamics 
of the magnetic moment in nanomagnetic materials at low temperatures, it giving rise to the ob- 
servation of quantum relaxations and quantum hysteresis cycles. For the first time, we report 
experimental evidence of quantum dynamics of the vortex core of micron-size permalloy (FeigNigi) 
disks induced by the application of an in-plane magnetic field. It is attributed to the quantum 
tunneling of the vortex core through pinning barriers, which are associated to structural defects in 
the dots, towards its equilibrium position. 

PACS numbers: 75.45.+j, 75.70.Kw, 75.78.Fg 
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Equilibrium magnetic configurations of soft ferromag- 
netic materials can be essentially non-uniform. It has 
been widely reported that, in the absence of applied mag- 
netic field, micron-size disks of these materials exhibit 
the vortex state (magnetic soliton) as the magnetic equi- 
librium of the system [lrQ|. Several technological and 
biomedical applications of the vortex state have been ex- 
plored, such as nonvolatile memory devices biomolec- 
ular carriers Q and targeted cancer-cell destruction 
The vortex state is characterized by a curling magneti- 
zation and an out-of-plane magnetic core, whose size is 
comparable to the material's exchange length (~ 6 nm). 
Striking feature of this magnetic equilibrium configura- 
tion is the unique role played by the vortex core: despite 
the fact it does not have any significant effect on the 
static properties of the system (it occupies far less than 
0.01% of the sample volume), the vortex core entirely 
governs the low frequency spin dynamics. In particu- 
lar, the excitation spectrum of these micron-size disks 
is characterized by the gyrotropic mode, corresponding 
to the spiral-like precessional motion of the vortex core 
as a whole 0, HI ■ Its frequency belongs to the sub-GHz 
range, and so it is intrinsically distinct from the spin wave 
modes. Furthermore, the chirality of this spiral motion 
is determined by the vortex core polarization p = ±1 

The vortex core is a suitable candidate to observe 
macroscopic quantum phenomena. Because of the strong 
exchange interaction it behaves as an independent entity 
and, the vortex core being a nanoscopic object, it is fea- 
sible that it exhibits quantum tunneling between clas- 
sically stable configurations. The measurement of time 
relaxations of the magnetic moment is a simple experi- 
mental way to observe this phenomenon. At finite tem- 
perature these relaxations may occur via thermal activa- 
tion, whereas in the limit T — > these relaxations con- 
tinue independently of the temperature due to underbar- 
rier quantum tunneling. Macroscopic quantum tunneling 



(MQT) i determines that the relaxation rate decreases 
as exp(— S e ff), with S e ff the action evaluated at the 
magnetic thermon (instanton), which takes into account 
the dynamics of the vortex core ruled by an energy po- 
tential and the dissipation of the system at a given tem- 
perature. This behaviour has been widely observed in a 
large number of systems [l(| , which include single domain 
particles fill. [Till , magnetic clusters [13], magnetic do- 
main walls |l4l|7fiux lines in type-II superconductors [HI 
and, very recently, Normal-Superconducting interfaces in 
type-I superconductors (lrj |. All these experimental evi- 
dences suggest that magnetic tunneling is a common phe- 
nomenon characterizing the low-temperature dynamics 
of magnetic materials in the mesoscopic scale. 

The application of an in-plane magnetic field yields 
the displacement of the vortex core perpendicularly to 
the field direction Q (see Fig. lb). It has been previ- 
ously reported that the dynamics of the vortex core can 
be affected by the presence of structural defects in the 
sample 17], |l8|. The gyrotropic motion of vortices, being 
the softest dynamical mode, is the most likely candidate 
for quantum relaxation. In the present letter we explore 
the magnetic irreversibility and the dynamics of vortex 
cores in micron-size permalloy dot arrays at low tempera- 
tures by means of the application of an in-plane magnetic 
field. For the first time, we report experimental evidence 
of the quantum depinning of magnetic vortex cores. 

All measurements were performed in a commercial 
SQUID magnetometer capable to measure at tempera- 
tures down to 2K and to apply magnetic fields up to 5T. 
The system is equipped with a Continuous Low Temper- 
ature Control (CLTC) and an Enhanced Thermometry 
Control (ETC) and it showed thermal stability better 
than 0.01 K at all times in any isothermal measurement. 
We have studied an array of permalloy disks with diam- 
eter 2R = 1.5 fim and thickness L — 95 nm. Its surface 
density is 0.15 dots//im 2 . Fig. la shows an AFM image 
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of this array. The array of permalloy disks was fabricated 
on a silicon wafer using optical lithography, and lift-off 
techniques: A single layer resist spin coating and highly 
directional electron-beam evaporation under UHV condi- 
tions were used to obtain circular dots with sharp edges. 
Optical lithography allows low cost patterning of a large 
area dot array in a single step. Consequently, identi- 
cal properties of magnetic material, such as grain size, 
distribution, and orientation, and film thickness may be 
obtained over the whole array. The magnetic film was 
deposited on a water-cooled substrate from a permalloy 
(FeigNigi) target. The growth ratio was of 1.5 A/s. 
The 95 nm Py film showed a switching field of about 4 
Oe and appeared to be pretty isotropic (in-plane). Fi- 
nally, the sample was prepared by stacking four 5x5 
mm 2 of these arrays with parallel sides and all magnetic 
measurements were performed using an in-plane config- 
uration for the applied magnetic field. The sample was 
studied in the range of temperatures 2 — 300 K and under 
applied magnetic fields up to 2 kOe. 




FIG. 1: a, AFM image of the array of permalloy disks studied 
in this paper, b, Spin field of the vortex state in one of the 
permalloy disk considered in a. The vortex core is displaced 
transversely to the direction of the applied field H. 
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FIG. 2: M(H) loops obtained at different temperatures (2, 50 
and 300 K) for the range of positive applied magnetic fields. 
The inset shows their zoom in the field range from H — — 50 
Oe to H = 50 Oe. 

Fig. 2 shows the M(H) hysteresis loop of the sample, 
in the range of positive applied magnetic fields, at differ- 
ent temperatures. For the negative range, the cycles are 



antisymmetric. The first magnetization curves have been 
omitted. Notice that these hysteresis loops correspond to 
the single domain (SD)o Vortex transitions 0,0]. As the 
temperature is lowered, the nucleation field H n decreases 
and the annihilation field H an increases (as reported in 
ref. [Hj]). For the range of temperatures explored, the 
vortex linear regime in the ascending branch should ex- 
tend from —300 Oe to 500 Oe at least. Inset of Fig. 
1 shows that the descending and ascending branches do 
not overlap at any temperature, and this extends over the 
whole linear regime. Furthermore, the remnant magneti- 
zation increases when T decreases [13] • In conclusion, the 
vortex linear regime exhibits magnetic irreversibility and 
it is temperature dependent. Consequently, we proceed 
to explore the metastability of vortices by means of i) 
ZFC-FC curves (Mzfc and Mpc) at different magnetic 
fields, and ii) isothermal measurements of the magnetiza- 
tion along the descending branch of the hysteresis cycle 
(Md es (H)), from the SD state, at different T. In both i) 
and ii) the values of T and H at which the magnetization 
has been measured were the same. 
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FIG. 3: a, Temperature dependence of Mzfc (300 Oe) and 
M FC (300 Oe) in the range 2-150 K. b, Plot of M des (300 Oe), 
together with Mzfc (300 Oe) and Mfc (300 Oe), in the range 
2 — 30 K. See text for details. 



The ZFC process consists of first performing minor cy- 
cles around H = Oe at T = 150 K (in order to get a 
zero magnetization state at H = Oe). Secondly, the 
sample was cooled down to T = 2 K without applied 
magnetic field and, thirdly, a desired magnetic field, H, 
was applied. Then, the ZFC magnetization curve was 
measured from 2 K to 150 K. Sweeping back the tem- 
perature to 2 K we follow the FC curve. Fig. 3a shows 
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the ZFC-FC curves obtained at H = 300 Oe. The mag- 
netization increases strongly from an initial value at 2 
K to a maximum in the ballpark of T ~ 30 K. Then 
it decreases smoothly and reaches a plateau. The de- 
pendence of the FC curve on T is similar to that of the 
ZFC case at high temperatures but with slightly higher 
values of M . In the ballpark of T ~ 20 K, the magne- 
tization of the FC rises strongly, reaching its maximum 
value at T — 2 K. This temperature dependence of both 
Mzfc(H) and Mpc(H) is characteristic of the range 
of applied magnetic fields in the linear vortex regime. 
Additionally, isothermal magnetic measurements along 
the descending branch of the hysteresis cycle (Md es {H)), 
from the SD state (H = 1 kOe), have been measured 
at the same T than the ZFC-FC curves. Fig. 3b shows 
Mdes(300 Oe) obtained in the range 2 — 30 K, together 
with M ZFC (300 Oe) and M FC (300 Oe). The values of 
Mdes (300 Oe) decrease strongly when T increases in the 
range 2 — 20 K and above T ~ 30 K tend smoothly to 
the FC curve (not shown). The divergence between the 
Mzfc, Mpc and M<i es curves in the range T — 2 — 20 K 
indicates the existence of a strong magnetic irreversibility 
in this region, and therefore we will focus on this range 
from now on. 

In order to confirm that the FC curve is the magnetic 
equilibria of the system, we performed two sets of mea- 
surements of the isothermal time evolution of the mag- 
netization, M(T,t), when sweeping the temperature in 
increments of 1 K per 30 minutes, i) from 15 K to 2 K 
and ii) from 2 K to 15 K. The initial magnetic state for 
each set was prepared by means of the above ZFC process 
to the desired temperature, followed by the application 
of a magnetic field H = 300 Oe. In i) it is only observed 
magnetic relaxation of the sample at 15 and 14 K, which 
quickly reaches a stable value corresponding to the FC 
one. From this point on, sweeping the temperature down 
to 2 K only leads to a variation of the magnetization of 
the sample following the values of the FC curve, ii) shows 
that, in the whole range of temperatures, the magnetiza- 
tion relaxes. The initial value of each relaxation follows 
the time evolution of the previous one. Moreover, the 
amount of relaxed magnetization is approximately the 
same for T = 2 — 9 K and it decreases progressively for 
10 — 15 K with magnetization values tending to the FC 
ones. 

We explored the metastability of the system more 
deeply by performing relaxation measurements in the 
vortex linear regime. The amount of magnetization avail- 
able for relaxation is M(0) — M eq , where M(0) is the 
initial magnetization value and M eq corresponds to the 
equilibrium magnetization. On account of this (see Fig. 
3b), we will focus our study on the relaxation measure- 
ments of vortices from the metastable states of the de- 
scending branch. Fig. 4a shows the normalized irre- 
versible magnetization (left term of equation (1)) vs. Inf 
curves measured for two different applied fields (H = 



and 300 Oe) at the same temperature (T = 2 K). Only 
below T ~ 15 K the magnetization of the sample fits very 
well a logarithmic time dependence. In this range of tem- 
peratures, the magnetic viscosity S(T) of the sample can 
be calculated by means of the theoretical formula (loj 
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Fig. 4b shows the viscosity, as a function of T, for two 
different magnetic fields. In both curves we observe a 
plateau at low temperatures and, what is more impor- 
tant, it does not extrapolate to zero in the limit T — > 0. 
In the ballpark of T ~ 7 K, the viscosity increases up to 
the temperature ~ 11 K, from which it decreases. Re- 
laxation measurements of the ZFC state with different 
applied magnetic fields were also performed, obtaining 
similar results for the viscosity. 



^1,00 
^ 0,98 
§- 0,96 

>: 0,94 

^"0,92 
<r- 0,90 
^ 0,88 



■ H=0Oe 
• #=300Oe 




5,0 5,5 6,0 6,5 7,0 7,5 8,0 

ln(/) 

0,04 



0,03 
CO 0,02 
0,01 



0,00 



-i — . — i — i — i — i — i — . — i — i — i — i — i — . — i — i — i — . — |- 



•iii 



• /V= 300 Oe 
■ H=0Oe 



i 



_i i i i i i i i i i i i i_ 



2 4 6 



10 12 14 16 18 20 
7-(K) 



FIG. 4: a, Normalized irreversible magnetization vs. Int 
curves measured at T = 2 K for H = Oe and 300 Oc. 
b, Temperature dependence of the magnetic viscosity S(T) 
at H = Oe and 300 Oe. See text for details. 

A logarithmic time dependence of the magnetization 
in relaxation measurements indicates the existence of a 
broad distribution of energy barriers U in our system. 
Classically, these energy barriers can be overcome by 
thermal activation, whose probability is proportional to 
the Arrhenius factor exp(— U /T). The so-called Blocking 
temperature, Tb, sets apart both reversible (T > Tb) 
and irreversible (T < Tb) regimes when the sample is 
externally perturbed. Despite the slight differences be- 
tween Mzfc and Mpc and between Md es and Mpc at 
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high temperatures, the strong divergence of magnetiza- 
tion observed in these curves suggests that the Blocking 
temperature should be below T ~ 20 K. 

Conventionally, the Blocking temperature for weakly 
interacting systems can be estimated as the tempera- 
ture at which the magnetic viscosity reaches its maxi- 
mum [10| . From our data we estimate that Tb ~ 11 
K, which is in good agreement with the gradual loss of 
logarithmic time dependence of our relaxation measure- 
ments at T > 15 K. Notice that thermal activation of 
energy barriers dies out in the limit T —¥ 0. Therefore, 
our observation that magnetic viscosity S(T) tends to a 
finite value different from zero as T — > indicates that re- 
laxations are non-thermal in this regime, i.e., transitions 
from metastable states are due to underbarrier quantum 
tunneling. This interpretation is upheld by the M vs. 
Tln(£/To) graphic representation [12j. The time r is the 
so called characteristic time attempt of the system and 
we have estimated its value to be To ~ 1CP 11 s, so that 
all magnetic relaxation curves only scale for temperatures 
above T ~ 9 K (see Fig. 5). This loss of scaling corre- 
sponds to the quantum regime case and is independent 
of the energy barrier distribution. 
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FIG. 5: Magnetization vs. Tint curve measured at H — 300 
Oe. Above T ~ 9 K, we verify the scaling M = M(T ln(f/r )), 
which corresponds to the case of thermal relaxation. Below 
T ~ 9 K we find a breakdown of this scaling, which corre- 
sponds to the quantum regime case. 



The onset of irreversibility appears sweeping the ex- 
ternal magnetic field for both the ZFC curve and the 
descending branch. The effect of this field is to move 
the vortex core along the disk surface and the observed 
irreversibility should come from this movement. Recent 
experimental data report the existence of some sort of 
structural defects in the disks [ijj , which could be a feasi- 
ble origin of the energy barriers responsible for the mag- 
netic dynamics of the system. In the light of this, we 
consider these defects are capable of pinning the vortex 
core. Therefore, the relaxation of the sample could be 
understood as simply the dynamics of the vortex core 



when escaping from the pinning centers and overcoming 
the corresponding energy barriers towards the equilib- 
rium. Note that lower temperatures are probing lower 
energy barriers because of the exp(— U /T) dependence of 
the probability of thermal activation. 

Within the framework of the rigid model of the shifted 
vortex, the vortex core is described as a zero-dimensional 
object whose dynamics is ruled by Thiele's equation. The 
corresponding Langrangian is given by C = Gyx — W(r), 
where r = (x, y) are the coordinates of the vortex core 
in the XY plane, G is the modulus of its gyrovector and 
W(r) is the total magnetic energy of the system. To in- 
corporate the smallest pinning barriers into the model 
we treat the vortex core as a stack of pancake vortices, 
one in each atomic layer. This pancake structure shows 
a finite rigidity in the vertical dimension, which means 
that these layers interact elastically among them. We 
consider that just a small vertical segment of the vortex 
core (equivalently, a small amount of these vortex layers) 
takes part in the tunneling of the vortex core through 
the pinning barrier, whose length is I << L. In order 
to model all this we consider the vortex core as a flex- 
ible line that goes predominantly along the z direction, 
so that r = r(z, t) is a field depending on the vertical 
coordinate of the vortex core, z. The whole magnetic en- 
ergy (including the elastic and the pinning potential) is 
described via a biparametric quartic potential given by 
W(r) = -^hx+\n{x 2 + y 2 ) - ±(3x 4 + f (ff)^, where 
/i and h are the magnetic moment of the dot, respectively 
the modulus of external magnetic field (which is applied 
in the y direction), A is the elastic coefficient of the pan- 
cake structure and k and (3 are the parameters of the 
potential energy. Assuming a second order transition for 
thermal to quantum relaxation, the obtained expressions 
for the crossover temperature T c and for the height of the 
barrier U in absence of applied magnetic field (h = 0) are 
the following [H 
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where the modulus of the gyrovector is given by the for- 
mula G = 2ir(+l)lM s /j, so that it is related to the tun- 
neling vertical segment. 

Comparison of the theoretical model with the experi- 
mental results leads to the determination of the param- 
eters (k,/3): first of all, notice that I cannot be smaller 
than the material's exchange length because, otherwise, 
the deformation of the vortex core line would be energet- 
ically unfavorable to the system. The same happens if I 
is much bigger than this exchange length. So it should 
be l~l ex = y / 2A/^ M 2 ~ 6 nm, where A = 1.3 • 1CT 11 
J/m is the exchange constant and M s — 7.5 • 10 5 A/m 
is the saturation magnetization of permalloy. The value 
of the modulus of the gyrovector is G = 2n(+l)lM s /j = 
1.62 ■ lCT 13 Ns/m (7 = 1.76 ■ 10 11 (Ts) _1 ). Experimen- 
tally, we have T c ~ 9 K for the H = Oe case too, from 
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which we deduce the value k ~ 0.5 J/m . On the other 
hand, for a measurable tunneling rate S e f / should be in 
the ballpark of 30. As S e jf — cnG/2^/2h(3 with c being 
a numerical factor of order unity resulting from the inte- 
gration, we have the following estimate of the coefficient 
P, P ~ nG/m^h = 9.8 • 10 18 J/m 4 . Finally, from these 
values of the parameters of the pinning potential we can 
estimate the width of the energy barrier, which is given 
by the expression Lb = y/2n/P ~ 0.3 nm. Furthermore, 
we can also estimate the order of magnitude of the height 
of the barrier (mean value), U ~ 250 K, which is in good 
agreement with the value given by the Blocking tempera- 
ture. These estimates are feasible values because pinning 
happens at the atomic level. 

In conclusion, the non-thermal dynamics of magnetic 
vortices in micron-size permalloy disks is reported. It 
is attributed to the quantum depinning of vortex cores 
through the structural defects of the sample, in steps 
about 0.3 nm. 
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